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A  -  SOURCE  DEPTH  DETERMINATION  AT  LOCAL  AND  REGIONAL  DISTANCES 


In  a  previous  study  (  see  Annual  Scientific  Report,  year  1982  )  we  tried 
to  evaluate  the  influence  of  source  depth  on  regional  phases  amplitudes 
recorded  between  20  and  500  km  of  epicentral  distances.  For  that 
purpose,  synthetics  were  computed  by  the  Bouchon  method  and  Pg  and 
Lg  waves  amplitudes  were  evaluated  for  various  depths  :  5,  17,  29, 

3  1  km. 

Besides  attenuation  1  s  evaluation,  it  clearly  appeared  that  Pg  waves 
amplitudes  are  rather  independant  of  depth  and  that  Lg  waves  amplitudes 
are  decreasing  versus  depth. 

Our  purpose  this  year  is  to  evaluate  if  this  method  could  help  to  discrimi¬ 
nate  events  versus  depth.  We  have  started  with  local  epicentral  distances 
(  1  0  to  60  km  ). 


Method  of  synthetics  computation 

The  method  of  synthetics  is  based  on  an  exact  discrete  wave  numbers  compu¬ 
tation  (Bouchon  1982  ).  We  briefly  recall  that  it  allows  to  evaluate  the 
total  field  representing  the  ground  motion  S  (w)  as  : 


3  3  _ 

S  (<*»)=  X  Z  Mij(w)  Xj(w)  lj 

i=l  j=l 


in  the  frequency  domain 


2. 


with  : 

-  ( 4 )  i  =  1,3  base 

-  Mjj  3x3  seismic  moment  tensor  describing  the  slope  of 

the  source  and  the  time  history  of  rupture. 

This  moment  tensor  Mjj  is  a  product  of  two  factors  : 

•  s  (u)  describing  the  rupture  versus  time 

•  mij  scalare  factor  describing  the  source  geometry 

(  in  our  case  a  pure  strike  slope  ) 

One  representation  of  the  source  function  s  (ui)  is  given 
by  Brune  : 

„/  >  s(o) 

VI  +  (f/fc)2* 

-  X  j  (w)  is  representative  of  the  propagation  and  depends  on  the 

epicentral  distance,  and  the  source  depth. 

The  propagation  model  we  have  used  defines  a  crust  model  for  France 
(  Darpa  Annual  Report,  1984)  currently  used.  No  anelastic  attenuation  has 
been  taken  into  account  (Q  = 

Different  seismograms  have  been  computed  every  10km  of  epicentral  dis¬ 
tance  from  10  to  60  km,  for  different  depths  :  1,  3,  6,  9,  12  and  27km 

within  the  crust. 

Results 

Computed  impulse  responses  have  been  convolved  by  a  source  function  of 
Brune  type  with  a  corner  frequency  fc  =  3  Hz  and  an  asymptote  defined 
by  0<  =  2  (  Fig.  1  ). 

The  evolution  of  signal  with  depth  is  clearly  seen  :  very  complicated  with 
a  source  near  the  surface,  essentially  composed  of  guided  waves  in  the 
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VELOCITY  SYNTHETIC  SISMOORAMS  AND  CORK ES PONDING  DISPLACEMENT  FREQUENCY  SPECTRA 


Fig.  1-e  -  VELOCITY  SYNTHETIC  SISMOGRAMS  AND  CORRESPONDING  DISPLACEMENT  FREQUENCY  SPECTRA 
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first  layers,  it  becomes  simpler  and  simpler  versus  depth  :  only  direct 
P  and  S  waves  still  exist. 

In  the  frequency  domain,  this  effect  seems  to  be  clearer. 

We  have  computed  amplitude  spectra  for  : 

-  Pg  phases  (time  window  defined  between  V  =  6,  5  km/s  an  d  3,  8km/s  ) 

-  Sg  phases  (  "  "  "  "  V  =  3,  8km/ s  and  2,  5km/ s  ) 

-  Rg  waves  (  "  "  "  "  V  =  2,  5km/ s  and  0,  8km/s  ) 

The  computed  spectra  are  shown  on  Fig.  2  (a,  b,  c,  d,  e,  f)  for  distances 
rounding  from  10  to  60  km. 

Sg  and  particularly  Pg  phases  are  rather  less  influenced  by  depth  except  for 
very  short  distances.  The  Rg  phases  (  surface  waves  in  sedimentary  layers  ) 
have  their  max.  to  min.  spectral  ratio  ranging  between  100  and  1  000.  The 
effect  of  distance  is  clearly  seen. 

In  order  to  be  independant  of  this  distance  effect  rough  synthetics  are  decon¬ 
volved  from  geometrical  attenuation. 


inversion  of  the  source 

inversion  technique  is  applied  in  the  frequency  domain  to  the  synthetics  to 
obtain  simultaneously  the  source  spectrum  and  the  geometrical  attenuation 
this  being  done  on  each  phases  for  each  distance  and  each  depth  previously 
selected. 

Geometrical  attenuation  coefficients  are  largely  .ncreasing  with  depth  for 
Sg  waves  from  1,25  (1  km  )  to  2,  50  (  6  km  )  and  then  rather  independant 
(9  to  12  km  ). 

For  Pg  waves,  geometrical  attenuation  seems  slightly  dependant  on  depth, 
increasing  from  1 ,  70  (1  km )  to  2, 00  (6km)  and  decreasing  to  1,70 
(  1  2  km  ). 


SYNTHETIC  DISPLACEMENT  SPECTRA  AT  ONE  EPICENTRAL  DISTANCE  FOR  VARIOUS  SOURCE  DEPTHS 


EPICENTRHL  DISTANCE 


Fig.  2-b  -  SYNTHETIC  DISPLACEMENT  SPECTRA  AT  ONE  EPICENTRAL  DISTANCE  FOR  VARIOUS  SOURCE  DEPTHS 
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SYNTHETIC  DISPLACEMENT  SPECTRA  AT  ONE  EPICENTRAL  DISTANCE  FOR  VARIOUS  SOURCE  DEPTHS 
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SYNTHETIC  DISPLACEMENT  SPECTRA  AT  ONE  EPICENTRAL  DISTANCE  FOR  VARIOUS  SOURCE  DEPTHS 


EPICENTRRL  DISTANCE 
S0  km 


SYNTHETIC  DISPLACEMENT  SPECTRA  AT  ONE  EPICENTRAL  DISTANCE  FOR  VARIOUS  SOURCE  DEPTHS 


EPICENTRflL  DISTANCE 
60  km 


SYNTHETIC  DISPLACEMENT  SPECTRA  AT  ONE  EPICENTRAL  DISTANCE  FOR  VARIOUS  SOURCE  DEPTHS 
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For  Rg  waves,  geometrical  attenuation  is  less  effective  but  also  depending 
on  depth  :  1 , 40  (  1  km  )  increasing  to  1 , 90  (3  km  )  and  decreasing  to  0,  50 

(  12km).  At  larger  depths  these  waves  are  not  excited  and  the  geometrical 
attenuation  coefficient  has  not  more  signification. 

On  the  Fig.  3  are  given  the  computed  source  spectra  for  each  phase,  at  dif¬ 
ferent  depths  and  the  corresponding  Brune  spectra  obtained  by  least  squares. 
Large  differences  between  the  original  sources  (  Brune  model  )  and  the 
sources  computed  by  inversion  are  observed  especially  in  their  low  frequen¬ 
cy  part.  They  are  not  interpreted. 

Nevertheless,  if  we  represent  (  Fig.  4  )  amplitude  ratio  of  these  sources  two 
by  two  (S/R,  P/R,  P/S)  versus  depth  term  we  computed,  we  observe  a 
sensible  variation  of  these  three  numbers.  (Amplitude  is  the  low  frequency 
mean  value  of  the  spectrum). 

Although  ditferent  points  are  still  not  be  solved  as  the  absence  of  flat  fre¬ 
quency  response  of  each  computed  source  spectrum,  we  might  expect  to 
determine  source  depth  ranges  at  local  distances  by  a  careful  analysis  of 
the  relative  importance  of  eaci'  nhase. 


Conclusion 

By  using  synthetics  computed  for  different  source  depths  and  for  a  distance 
range  from  10  to  60km,  the  source  function  has  been  inversed  for  different 
phases  (  Pg,  Sg  and  Rg  phases)  defined  by  their  group  velocity;  the  pro¬ 
pagation  being  caracterized  by  its  geometrical  attenuation  only. 

Besides  a  larj.c  dispersion  of  the  results  and  the  fact  that  it  is  not  always 
possible  retrieve  the  original  source  function  used  to  compute  the  synthetics, 
it  seems  that  the  low  frequency  part  of  the  displacement  source  spectra  we 
obtained  is  strongly  influenced  by  depth. 
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SOURCE  SPECTRA  ( Depth-  1  Km) 


II:  DOTTKI)  r  ,IME  ( lot;  lot;  scale)  : 

corner  frequency  -  <  *1 1  111*’  slope  -  A  ;  tioometrical  attenuation  factoi 


SOURCE  SPECTRR  (Depth-  3  Km) 


SOURCE  SPECTRR  (Depth-  6  Km) 


BONDING  BRUNE  SPECTRA  IN  DOTTED  LINE  (log-log  scale)  : 

frequency  amplitude  -  fc  =  corner  frequency  -  ‘X  J  HF  slope  -  A  :  geometrical  attenuation  factor 


SOURCE  SPECTRR  (Depth-  9  Km) 


Fig.  i-e  -  CORRESPONDING  HKIJNE  SPECTRA  IN  DOTTED  LINE  (log-log  scale)  : 

Ao  :  low  frequency  amplitude  -  fc  :  corner  frequency  -  *  \  :  II  F  slope  -  A  :  geometrical  attenuation  factor 


,OW  FREQUENCY  SPECTRAL  RATIO  VERSUS  DEPTH 
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If  that  is  confirmed,  we  might  obtained  information  on  the  depth  of  a 
seismic  event  at  least  at  local  distances  by  using  the  ratio  of  Rg  upon 
Pg  low  frequency  spectra. 

Nevertheless,  it  is  necessary  to  verify  that  the  curve  we  obtained  is 
independant  on  the  other  source  characteristics  that  is  the  source 
mechanism. 


B  -  INVERSION  METHOD  ON  Q  ( f )  -  Q  TOMOGRAPHY  ON  FRANCE  AREA 


The  studies  we  are  carrying  on  to  interprete  the  regional  seismic  phases 
in  term  of  source  energy  have  pointed  out  the  main  effect  of  anelastic 
attenuation. 

After  having  studied  the  propagation  of  regional  phases  on  a  theoretical 
basis  by  synthe,.  .3,  we  have  worked  on  the  quantitative  interpretation  of 
the  seismograms. 

The  data  collected  from  the  L.  D.  G.  network  have  led,  in  a  first  step,  to 
represent  the  anelastic  attenuation  by  a  mean  quality  factor  depending  on 
frequency. 

This  model  has  been  confirmed  by  the  coherency  between  empirical  data 
and  numerical  simulation. 

Consequently,  it  was  possible  to  compute  the  source  excitation  function 
from  each  seismogram.  The  residual  dispersion  of  source  functions 
obtained  by  this  way  might  be  explained  by  a  spatial  heterogeneity  of  the 
quality  factor  in  the  France  region. 

This  heterogeneity  is  associated  to  different  geological  units  with  proper 
quality  factor  and  typical  frequency  dependancy. 

As  a  result,  we  have  used  tomography  techniques  to  evaluate  the  local 
quality  factor. 

A  "back  projection"  application  has  pointed  out  a  strong  attenuation  hete¬ 
rogeneity  at  1  Hz  which  geographical  di stnbution  is  strongly  correlated  with 
local  geological  characteristics. 

These  different  results  have  led  us  to  try  to  improve  our  seismic  energy 
measurements  from  Lg  waves  analysis  by  taking  into  account  the  propagation 
anomalies. 
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For  that  purpose,  it  is  necessary  to  draw  a  detailled  map  of  the  quality 
factor  by  an  inversion  technique  at  each  frequency 

Such  an  operation  needs  a  large  number  of  data  which  are  only  available 
from  natural  sources. 

The  source  energy  at  each  frequency  is  becoming  an  extra -unknown  for 
which  the  specific  response  of  each  station  is  to  be  added. 

This  study  will  allow  an  evaluation  of  the  real  resolution  capacity  in  ener¬ 
gy  given  by  an  homogeneous  set  of  regional  phases  records. 

The  large  amount  of  data  and  unknowns  leads  to  use  a  technique  of  iterative 
reconstruction  which  needs  not  the  resolution  of  large  equation  systems. 

We  have  developped  a  program  of  local  reconstruction  (SIRT)*.  This 
technique  appeared  to  be  well  adapted  to  solve  the  very  inhomogeneous 
spatial  distribution  between  the  source  and  each  station. 

We  have  introduced  in  the  reconstruction  algorithm  a  weighting  function  taking 
into  account  the  possible  redondancy  of  information  due  to  the  peculiar  geo¬ 
graphical  distribution  of  natural  seismicity.  The  weighting  function  applied 
for  each  seismic  path,  at  each  point  is  depending  on  path  azimuth. 

We  have  tested  our  program  on  synthetic  models  and  on  limited  data  sets  and 
are  just  to  proceed  on  a  routine  basis  on  a  large  amount  of  data  recorded 
by  the  L.  D.  G.  network  and  as  we  expect  by  other  european  digital  networks 
(Switzerland,  Italy). 


Simultaneous  Iterative  Reconstruction  Technique 
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C  -  SOURCE  STUDIES  -  DATA  BASE 


In  1986,  we  have  achieved  a  regional  seismic  data  base  which  includes 
the  seismic  data  of  431  european  earthquakes  recorded  on  the  french 
L.D.G.  network  (Fig.  1),  during  1984  to  1986,  with  ML  magnitudes 
ranging  from  2.  4  to  5.  7  (  Fig.  2  ).  All  these  quakes  are  within  the 
crust  but  for  most  of  them  depth  is  not  constrained. 

The  data  base  contains  the  seismic  records  (samples  at  50  samples/sec) 
and  the  source  spectra  of  each  quakes  computed  on  Lgj  Lg2  Lgj  and 
Pg  waves  inversion  techniques  (see  annual  scientific  report  -  1985). 

In  order  to  select  records  with  clear  Lg  waves,  seismic  signals  are  re¬ 
read  from  numerical  mag-taps  on  paper  records  with  a  zero  phase  on 
Pg  waves  (  Fig.  3  as  an  example  ). 


Lg  waves  train  is  divided  in  three  parts  according  to  the  group  velocity  : 


3 ,  1  km/  s  < 

V  < 

3 ,  6  km/  s 

2 ,  6  km/  s  < 

V  < 

3 ,  1  km/  s 

2,  3  km/  s  < 

V  < 

2,  6  km/  s 

For  each  record,  the  amplitude  spectrum  is  computed  for  each  of  the  three 
Lg  phases  between  0,  5  and  15  Hz.  The  seismic  noise  is  also  computed 
in  front  of  each  signal  integration  time  in  order  to  appreciate  the  usefull 
frequency  oand  for  the  three  Lg  phases. 

The  whole  set  of  spectra  is  the  data  base. 

In  order  to  test  the  data  validity  a  first  study  on  source  spectra  has  been 
conducted  by  deconvolution  of  Lg  spectra  from  the  homogeneous  attenua¬ 
tion  model  we  obtain  previously  (  Campillo  and  al.  ,  1985). 


Fig. 


NUMBER  OF  EARTHQUAKES  VERSUS  ML  MAGNITUDE  HISTOGRAMS 
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For  each  quake,  the  source  spectrum  is  obtained  from  each  phase  Lg] 
Lg2  and  Lg 3  ■  and  also  for  the  combinde  phases  set  (  Fig.  4  ). 

A  large  amount  of  such  spectra  computed  in  each  station  have  a  large  dis¬ 
crepancy  which  might  be  reduced  in  a  reasonable  scale  by  the  introduction 
of  a  Q  tomography. 

When  the  dispersion  on  deconvolved  spectra  is  not  too  larger  -  which  is 
true  for  50  %  of  cases  -  their  mean  value  leads  to  obtain  the  source 
amplitude  at  a  distance  of  1  km  and  consequently  an  estimation  of  the  seis¬ 
mic  moment  computed  on  Lg  waves. 

The  source  expression  is  defined  as  a  Brune  source  of  the  form  : 


s  (w) 


«(o) 


\/l  +  (f/fc)2* 


and  the  sources  are  computed  by  least  squares  to  fit  this  amplitude  spectrum. 

For  the  different  parts  of  the  Lg  spectrum  the  seismic  moment  M0  is 
systematically  computed  and  the  relationships  between  M0  and  ML 
local  amplitude  are  built  up  (Fig.  6,  7,  8). 

W  e  obtain  : 


M°Lgl  » 

15,0 

+ 

1.4 

M°  Lg2  = 

15,  5 

+ 

1,3 

ml 

M°Lg3  = 

15,  6 

+ 

1,2 

ml 

For  the  whole  set  of  Lg  phase  : 

MLg  =  15.  1  +  1,4  Ml 


Lg  = 


10H2 


CM- 


i 

i 


OSSAU 

190486  22  5  13.6  43.06  -0.45  3.3 


33. 


1 


34. 


As  a  result,  the  seismic  moment  variation  versus  magnitude  is  very 
similar  for  each  part  of  the  spectrum,  according  to  the  large  discrepancy 
on  both  terms  aj  and  bj  from  : 

MLgj  =  ai  +  t>i  Ml 

This  relationship  MLg  =  f  (  Mj^  )  is  compared  with  the  relationships 
obtained  in  other  parts  of  the  world  with  close  in  epicentral  seismic 
data 


Mo  = 

16,2  + 

1,  37  Ml 

Japon 

M0  = 

17,  5  + 

1 , 20  Ml 

Californie 

M0  = 

16, 97  + 

1,  17  Ml 

Yunam 

For  France  area,  the  a  term  we  computed  is  1  or  2  units  weaker  than 
the  previous  ones,  this  might  be  interpreted  as  a  source  depth  effect  which 
was  not  be  taken  into  account  in  our  computation. 

An  attempt  to  regionalize  quakes  by  peculiar  M0  =  f  (Ml)  relationship 
did  not  lead  to  positive  results.  Q  tomography  might  bring  more  accurate 
attenuation  terms  and  so,  allow  some  quakes  regionalization  if  any. 
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In  the  last  annua!  report  we  have  described  a  method  (  Campillo  and  ai  1984) 
to  extract  from  the  Lg  waves  train  the  Q  factor  associated  to  the  propagation 
region  and  the  source  spectrum. 

We  have  applied  these  results  for  some  quakes  occuring  in  Remiremont 
(  Vosges)  in  December  84/ January  85  and  recorded  on  the  seismic  stations 
of  the  LOG  network.  The  aim  of  this  study  was  to  evaluate  the  source 
amplitude  spectrum  for  the  main  shock  and  some  precursors  anc  after 
shocks,  and  check  if  these  results  extracted  from  Lg  waves  in  the  far 
field  were  coherent  with  those  obtained  at  short  distance  from  the  epicenter. 

Seismic  data 

Among  the  28  SP  seismic  stations  of  the  LDG  network  only  8  having  recorded 
this  crisis  with  a  reasonable  signal  upon  noise  ratio  have  been  selected. 

In  the  Figure  1  are  represented  the  data  recorded  for  the  main  shock 
(  12^29/84  ;  Ml  -  4.8  )  and  the  stations  versus  epicenter  repartition. 

All  quakes  are  superficial,  between  5  and  10  km  depth.  Seismic  data  are 
numerically  recorded  on  magnetic  tape  with  a  sampling  rate  of  the 
50  samples  /  sec  which  allows  to  process  signals  spectra  up  to  15  Hz. 
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Results  : 

-  After  shock  01,02.85  ML  =  4.0 

This  after  shock  is  recorded  at  HAU  station  (  ^  =  21  km  ) 

without  saturation. 

After  deconvolution  of  the  seismograph  response,  the  amplitude  spectrum 
of  the  Sg  signal  is  computed. 

It  represents  an  amplitude  spectrum  of  the  source. 

A  source  displacement  model  of  the  form  : 


n  (o 


’.9  computed  by  leas:  squares  to  fit  this  amplitude  spectrum  (  Fig.  2  a  ) 
We  find  :  corner  frequency  fc  =3.6  Ha 

slope  if  1 

The  seismic  moment  Mo  is  also  available  by  the  following  relation  ship 
(  Street  and  al  1975  ;  Hermann  and  Kijko  1983) 


n(o) 

f> 

r  : 
r6* 

We  find 


w  e  ft  3  ron(o) 

R  84> 

km 

displacement  source  spectrum  at  very  low  frequency 
shear  wave  velocity  =3.5  km/s 
density  2.2  g/cm3 

coefficient  taking  into  account  both  radiation  pattern  #  1 
Mo  =  2.8jtl0  ^  dyne  .  cm 


We  also  compute  from  Lg,  .  Lg  2  .  Lg  3  waves  trains  recorded 
in  the  8  stations  situated  at  more  than  200  km  from  the  epicentral  area, 
a  mean  amplitude  source  spectrum.  For  this  computation  we  correct 
the  Lgj  f  Lg  2.  Lg3  ,  spectra  from  : 

-  elastic  and  anelastic  attenuation 

-  station  response 

-  seismograph  response 


45. 


Figure*  2b,  2  c  ,  2  d  represent  an  evaluation  of  the  +  limit* 

for  the  8  eource  apectra  and  the  mean  amplitude  aource  spectrum. 

(  continue  line  ). 

Wt  also  figure  the  aource  spectrum  computed  from  HAU  Sg  waves 
(  doted  line  )  as  a  reference. 

Asa  result  ,  source  spectra  obtained  with  Lg  waves  in  the  far  field  are 
similar  to  the  one  obtained  in  the  near  field  (  HAU  )  namely.  Lg]  waves 
lead  to  a  source  spectrum  which  fl  (o),  fc  and  slope  are  very  comparable 
to  those  of  HAU. 

Seismic  moments  computed  on  the  source  spectra  from  Lg  waves  and  the 
one  extracted  from  Sg  waves  at  HAU  are  of  the  same  order. 


-  Main  shock  12.  29, 84  Ml  *  4. 8 

Source  amplitude  spectra  on  Lg  1  ,  Lgj  ,  Lg  3  waves  are  computed 
by  the  same  method  described  in  the  first  case  (  After  shock  01  02  85 
Ml  *  4.0). 

For  this  main  shock  case  no  records  are  available  at  HAU. 

Figure  3  represents  the  source  spectra  we  obtain  : 

-  -  O-  amplitude  spectra  corresponding  to  all  the  8  stations 

the  source  spectrum  which  fits  the  mean  amplitude  spectrum  values, 
corner  frequency  fc  ,  slope  of  the  spectrum  and  seismic  moments 
which  are  very  similar  . 

-  fault  dimensions  and  stress  drops  computed  for  each  of  the  three  cases 
with,  for  the  circular  fault  radius, two  different  formula  proposed  by  Brune 
and  Madariaga. 


L 


0.  37 


(  Brune  ) 


L  r  0 .  2 1  _ 

fc 

(  Madariaga  ) 


A  _ 
Ue-  = 


Mo 

~r 


and 
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-  Small  foreshock  12,24.84  Ml  *  3.0 

This  last  example  of  source  evaluation  it  given  to  figure  out  the 

limits  of  our  far  field  method.  This  small  foreshock  is  recorded  by  8  stations 

among  those  which  recorded  the  previous  cases. 

Figure  4  presents  the  source  amplitude  spectra  obtained  at  HAU  on 
Sg  waves  and  on  Lg]  ,  Lg2  ,  Lg3  waves  recorded  on  the  8  stations 
at  200  km  and  more  from  the  epicenter.  Computations  and  representations 
are  identical  as  the  ones  of  Figure  2  on  the  aftershock. 

For  this  small  foreshock  the  corner  frequency  is  of  the  order  of  7  to  9  Ha, 
to  be  compared  to  3  .  6  Hz  for  the  after  shock  (  Ml  =4.0)  and  around 
2  Hz  for  the  main  shock  (Ml  =  4.8  ). 

19 

The  seismic  moment  is  of  the  order  of  some  10  dyne  .cm. 

But  essentialy  ,  the  spectra  definition  particularly  for  Lg3  waves  and  even 
Lg2  is  here  rather  poor  and  do  not  guarantee  a  reasonable  fiability  on  the 
results . 


Conclusion  : 

This  attempt  to  evaluate  the  source  amplitude  spectrum  of  the  earthquake 
by  processing  the  Lg  waves  train  recorded  at  large  distances  seems  to  be 
reliable  as  shown  on  the  Remiremont  se;smic  crisis. 

It  is  necessary  to  apply  this  method  on  more  different  cases  but  we 
notice  already  that  as  long  as  the  Lg  waves  are  well  defined,  they  bring 
information  from  the  source  which  lead  to  approach  its  main  parameters. 
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Introduction 

The  method  used  in  that  study  is  derived  from  the  discret  waves 
numbers  method.lt  has  been  used  to  characterize  the  effects  of 
different  large  scale  heterogeneities  on  the  Lg  waves  propagation. 

Two  types  of  effect  have  been  studied. 

First,  the  effects  due  to  Moho  irregularities.  The  Lg  waves  tram 
is  a  superposition  of  reflected  S  waves  trapped  within  the  crust. 
Consequently  any  variation  of  the  Moho  depth  might  affect  the  waves 
guide  quality  and  leads  to  some  apparent  strong  attenuation  zone  for 
Lg  waves. 

Second,  we  look  at  the  effects  of  sedimentary  basins  on  the  Lg  waves 
propagation  path. 

The  study  is  limited  to  the  SH  in  two  dimensions.  Irregularities 
dimensions  are  selected  to  fit  the  situations  met  in  central  part  of 
France  for  which  general  attenuation  features  have  been  studied 
(  Campillo  ,  Planter  and  Bouchon  1985). 
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I  Variation  of  the  Moho  depth 

The  crust  model  is  given  in  Figure  1.  Computations  have  been  made 
for  seismic  sources  at  10  km  depth,  located  at  0,  50,  100,  150  and  200  km 
from  the  Moho  irregularity. 

Figure  2  shows  maximum  amplitude  variations  for  Lg  waves  versus 
distance,  for  a  planar  crust  mode  (  thick  solid  line)  and  for  a  crust 
model  with  irregularity  at  these  previous  distances. 

In  the  irregularity  vicinities,  steepen  reflectors  bring  strong  amplitude 
variations . 

On  the  contrary,  at  larger  distances,  no  important  and  systematic 
ditlerences  are  seen  between  a  Moho  without  or  with  irregularity. 

This  points  out  that  any  Moho  ascent  should  commonly  bring  limited 
attenuation  effects  on  lg  waves. 

II  Sedimentary  basin 

We  look  now  at  the  influence  of  a  thicken  sedimentary  layer  or  a  sedi¬ 
mentary  basin  as  described  in  Figure  3  on  Lg  waves  at  the  same  distances 
from  the  irr  egula  r  ity  fa  s  describe  in  Figure  1  )  seismic  sources  are  always 
at  I  0  km  depth. 

The  synthetic  seismograms  are  presented  in  Figure  4  for  various 
distances  from  0  to  290  km,  every  10  km,  for  both  a  planar  crust  (  upper 
seismograms  )  and  a  thicken  sedimentary  layer  at  the  top  of  the  crust 
(  lower  seismograms). 

Clearly  seen  is  the  formation  of  a  Love  wave  in  the  case  of  a  sedimentary 
basin.  As  superficial  layers  are  propagation  structures  which  strongly 
affect  Love  waves,  the  converted  Lg  waves  to  Love  waves  energy  will  be 
rapidly  dissipated. 

Within  the  sedimentary  basin  itself,  classical  amplification  phenomena 
might  happen. 


These  two  rough  synthetic  examples  suggest  to  verify  their  conclusions 
for  some  real  cases. 


The  set  of  seismic  data  which  leads  to  the  evaluation  of  the  quality 
factor 

Q  =  Q  „  f 

for  the  center  of  France  in  a  recent  study  is  used  to  build  up  an  image 
of  the  Q  factor  spatial  distribution. 

The  method  we  use  is  the  ART  (  Algebric  Reconstruction  Technics). 

The  set  of  used  propagation  paths  for  the  1  Hz  reconstruction  (that  is 
to  say  Q0  )  is  presented  in  Figure  5. 

Each  point  of  the  image  is  taking  into  account  the  paths  within  a  distance 
of  50  km. 

In  order  to  apply  some  confidence  criterium  on  the  quality  of  the 
Q  evaluation  we  define  a  confidence  index  : 

ic  =  n  source  x  ZI  *r 

where  :  n  source  is  number  of  sources  which  spectra  are  different, 
used  for  the  Q  measurement. 

lr  is  the  ray  path  lenth  within  the  influence  circle  (  of  50  km 
radius  at  the  most  )  divided  by  the  influence  radius. 

A  map  of  these  factors  ic  is  shown  in  Figure  6.  All  the  represented 
points  fit  the  length  of  ray  test. 

All  the  blue  zones  correspond  to  a  high  confidence  index;  red  zones  might 
include  systematic  bias. 

Figure  7  represents  a  rather  simple  regional  distribution.  If  one  just  keeps 
within  high  confidence  limits,  three  main  regions  are  defined. 

-  Limagne  Basin,  Forez,  Saone  valley  weak  figure  for  Q 
(  between  180  and  240  ) 

-  Limousin  Q  between  200  and  260 

-  South  parisian  basin.  Berry,  Poitou,  North  of  Massif  Central 
Q  betw  een  300  and  440 

-  At  1  Hz  the  mean  value  for  Q  on  the  whole  studied  zone  is  2Q0  as 
mer.tionned  in  a  last  study. 


52. 


Conclusion 


These  results  pointed  out  the  importance  of  graben  zones  for  Lg  waves 
propagation.  Nevertheless,  at  this  stage  of  the  study  we  do  not  possess 
any  separation  tool  between  two  types  of  phenomena  which  could  create 
attenuation  : 

-  small  scale  diffraction  (  or  scattering  ) 

-  large  scale  geometrical  effect. 

A  study  of  frequency  dependancy  for  Q  should  help  in  the  interpretation 
of  our  data . 

On  the  other  hand,  developments  of  numerical  simulation  methods 
encourage  us  to  investigate  more  complex  zones  in  order  to  interprete 
the  quasi-extinction  phenomena  of  Lg  waves  we  have  observed  along  peculiar 
paths. 
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Figure  1  :  crust  model  with  Moho  irregularity  -  the  star  is  the 

seismic  source  -  seismograms  are  built  at  various  distances. 

Figure  2  maximum  amplitudes  of  Lg  waves  versus  distance  (  2D  ). 

Figure  3  crust  model  with  a  sedimentary  basin 

Figure  4  :  synthetics  obtained  with  a  planar  crust  model  (  high)  and  with 

a  sedimentary  irregularity  (  low  ) 

Figure  5  Rays  paths  used  to  build  up  a  map  of  Q  0  factor  for  France 

Figure  6  Map  of  confidence  index 

Figure  7  :  Map  of  France  (  Center  )  for  Q  factor  at  1  Hz 
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In  despite  of  the  title  of  this  presentation,  we  will 
not  present  a  study  of  source  functions,  but  a  way  to  obtain 
a  signal  averaged  over  the  seismic  stations  of  the  French 
Seismic  Network,  which  will  enable  us  to  study  the  source  in 
a  future  work. 


I .  Obtention  of  the  averaged  signal. 

A  classical  way  of  averaging  signals  over  a  seismic 
network  is  to  use  a  "Delay  and  Sum".  In  this  procedure,  a 
simple  stack  of  the  various  signals  are  done,  without  taking 
into  account  the  change  in  the  records  of  each  station. To 
overcome  this  problem,  we  decide  to  define  what  we  will  call 
the  best  averaged  signal  in  a  least  square  sens. 

The  technic  used  is  far  from  beeing  new.  Ke  suppose 
each  of  the  records  to  consist  of: 

-  a  common  signal  but  multiplied  by  a  station's  factor. 

-  a  noise  (mostly  wave  generated)  completely  decorrelated 

from  the  common  signal. 

The  signals  are  first  suppossed  to  be  delayed,  so  to  be 
in  phase,  but  we  will  come  back  later  to  this  point.  We 
compute  the  correlation  matrix  of  all  signals.  It  can  be 
shown  that  the  station's  factors  is  one  of  the  eigen  vectors 
of  this  matrix,  the  one  corresponding  to  the  largest  of  all 
eigenvalues.  This  introduce  an  undetermination:  an  eigen 
vector  is  allways  obtained  up  to  a  constant,  but  this  is 
obvious  on  physical  basis,  it  is  possible  to  multiply  the 
station's  factors  by  an  arbitrary  number,  if  in  the  same 
time,  the  common  signal  is  divided  by  the  same  value. 

The  common  signal  is  then,  the  sum  of  the  records, 
averaged  by  the  just  computed  station's  factors.  The  wave 
generated  noise,  or  residual  signal  in  a  given  station,  is 


obtained  by  a  difference  between  the  record  in  the  station 
and  the  common  signal. 

This  decomposition  of  the  signals  is  quite  similar  to 
the  covariance  analysis,  well-known  in  statistical  studies. 
The  main  difference  is  that  in  statistics,  the  decomposition 
of  the  signals  is  usually  done  on  a  larger  set  of  eigen 
vectors  than  just  the  first  one.  But,  in  our  case,  only  the 
first  as  a  meaning,  corresponding  to  a  signal  beneath  France 
and  representing  the  source  function  filtered  by  the 
propagation.  A  decomposition  on  a  larger  set  will  have  no 
physical  meaning,  the  coda  in  each  station,  are  only 
statisticaly  similar.  Of  course,  such  a  decomposition  has  to 
be  done  over  a  seimic  network  which  aperture  is  small  enough 
compared  tr  the  epicentral  distance. 

This  procedure  is  relatively  sensitive  to  the 
determination  of  the  arrival  time  in  each  station.  We  were 
greatly  helped  in  this  by  the  origin  of  the  events  we 
studied.  It  consist  in  a  "seismic  swarm",  a  serie  of  very 
similar  events  from  the  same  area  in  Eastern  Kazakh,  not  far 
from  Shagan  River.  The  signals  recorded  in  one  station  for 
the  various  events  are  more  similar  than  the  signals 
recorded  in  the  stations  of  the  French  Seismic  Network  for 
the  same  event.  So,  in  a  first  step,  we  computed  the  common 
signal  in  each  single  stations  for  a  set  of  10  events.  It 
was  easier  to  measure  the  arrival  time  on  this  averaged 
signal,  then  to  obtain  a  precise  arrival  time  for  each 
record  by  a  correlation  technic.  We  also  checked,  by 
iteration,  that  it  was  not  possible  to  improve  the  time 
determination. 

We  vary  the  length  of  the  time  window  used  in  the 
computation.  The  samplig  rate  of  our  records  is  50Hz,  with 
an  anti-aliasing  filter  at  16Hz,  in  fact  all  of  our  data 
were  low-pass  filtered  at  8Hz.  The  various  time  windows  were 
between  1024  points  (atout  20  seconds)  and  128  points  (2.5 
seconds)  without  real  change  in  the  station's  factors. 


II.  Magnitude  determination. 

The  first  interpretation  is  to  consider  the  common 
signal,  in  the  various  stations,  as  representing  the  signal 
arriving  beneath  France  and  filtered  by  an  average  crust, 
the  station's  factors  as  a  measure  of  the  sensitivity  of 
each  station,  and  the  residual  signal  as  athe  coda  (in  fact 
we  obtain  more  than  the  coda,  we  also  have  energy  in  the 
beginning  of  the  signal,  see  later). 


We  will  start  by  magnitude  considerations.  We  studied 
two  different  ways  to  measure  the  magnitude.  In  the  first, 
for  each  station,  we  computed  for  all  events  a  common  signal 
and  ,in  that  case,  the  event's  factors.  The  magnitude  beeing 
the  logarithm  of  the  wave  amplitude,  plus  a  distance 
correction,  and  as  all  events  are  practicaly  at  the  same 
distance,  the  logarithm  of  the  event's  factors  should  be, 
within  a  constant  correction,  the  event's  magnitude. 

Using  the  magnitude  determination  of  the  French 
Network,  we  found  the  best  relation  between  magnitude  and 
log  of  event's  factor  to  be  of  the  form: 

mb  ■=  a*log(EF)  +  b 

but  with  a  "a"  factor  not  equal  to  1.  We  give  in  annex  the 
list  of  all  the  results,  the  mean  values  are: 

a  -  0.837  +/- . 086 

b  -  5.952  +/-.035 


We  also  used  only  nine  of  the  events  to  recompute  a  and 
b  ,  to  calculate  the  magnitude  of  the  last  one.  Volontarily, 
we  select  the  nine  smallest  events  to  obtain  the  magnitude 
of  the  largest.  The  result  were: 

m  -  6.17  +/--03 

which  is  exactly  the  value  found  by  the  French  Network. 

We  intent  to  apply  the  same  technic,  on  a  much  larger 
data  set,  in  order  to  obtain,  not  really  a  magnitude,  but  a 
classification  of  the  various  events,  according  to  their 
energy. 

This  first  way  to  determine  a  magnitude  is  rather 
tedious:  each  time  a  new  event  is  aded  to  the  data  base,  all 
computations  have  to  be  redone.  An  other  way,  is  to  use  the 
mean  station's  factors.  We  noticed  that  the  station'  factors 
are  very  similar  for  a  given  station  and  the  various  events. 
So,  we  first  computed  an  average  value  of  the  station's 
factor,  in  each  station  for  all  events,  then  we  computed  an 
approximation  of  the  common  signal,  using  this  new  set  of 
station's  factors  and  measured  it's  amplitude.  The  results 
are  given  in  annex. 


III. 


Source  function. 


The  common  signal,  from  all  stations  of  the  network  and 
a  given  event,  represents  an  approximation  of  the  incoming 
signal,  beneath  France  and  filtered  by  the  mean  french 
crust.  We  present  here  some  of  the  common  signal  ve 
obtained.  In  this  preliminary  part  of  the  work,  we  have  not 
yet  compare  our  results  with  a  synthetic  sismogramm,  using  a 
source  function  and  a  propagation  effect. 

In  the  following  figures,  we  plotted  the  common  signals 
for  all  ten  events  used  in  this  study,  and  also,  for 
comparison,  the  actual  signals  in  some  stations  both  for 
the  best  event  and  an  average  one. 


IV.  Signal  residuals. 

The  difference  between  the  actual  record  and  the  common 
signal  represents  more  than  the  coda.  As  we  pointed  earlier, 
we  also  have  energy  in  the  first  part  of  the  signal. 

In  the  figures,  we  show  the  mean  signal  computed  using 
all  events  in  the  fives  stations  of  Massif  Central,  in  the 
middle  of  France.  And,  for  the  same  stations,  the  mean 
residual  signal. 

On  this  last  plot,  some  coherency  between  the  signals 
accross  the  subarray  is  clearly  seen. 

We  think  that  the  size  of  the  seismic  network  is  an 
important  parameter  for  the  coherency  of  the  waves  accross 
it.  If  the  aperture  is  too  small,  we  expect  to  obtain  a  high 
coherency  for  the  signals  themselves  (up  to  the  amplitude 
which  is  known  to  varry) ,  so  the  residual  signals  should 
represent  the  coda.  If  the  aperture  is  larger,  then  it  is 
possible  to  observe  propagations  through  the  array.  We  know, 
from  other  studies,  that  the  crust,  just  south  of  the  sub¬ 
array,  presents  larg  anomalies.  We  whish  to  correlate  both 
phenomens. 


V.  Preliminary  conclusions. 

This  presentation  is  only  a  presentation  of  a  current 
work.  Our  method  for  precise  magnitude  determination  has 
already  given  us  very  good  preliminary  results,  and  we  will 
improve  our  statistic  by  increasing  the  number  of  events. 
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The  common  signals,  obtained  for  a  given  quake,  over 
our  network  will  allow  a  study  of  the  source  function.  And 
we  really  hope  to  explain  the  main  trends  of  the  residual 
signals  by  correlation  with  the  geological  structure  around 
the  network. 


-JL. 


Annex 


Magnitude 


0.7527 

0.8390 

0.8184 

0.7185 

0.8635 

0.7973 

0.8387 

0.8313 

0.9131 

0.7698 

0.8024 

0.9291 

0.8525 

0.8399 

0.8257 

0.8985 

0.8849 

0.8236 

0.8575 

0.8890 

0.8103 

0.8578 

0.8799 

0.8114 

0.7131 

1.1354 

0.6380 


coef.  with  all  events 
b 

5.9658 

5.9690 

5.9445 

5.9974 

5.9124 

5.9924 

5.9021 

5.9924 

5.9124 

5.9687 

5.9529 

5.9509 

5.9552 

5.9574 

5.9629 

5.9497 

5.9441 

5.9285 

5.9212 

5.9227 

5.9494 

5.9210 

5.9135 

5.9365 

5.9693 

5.9337 

6.0732 


I 

t 


Magnitude  computed  from  common  signal 
m»a*log(A) +b 


a-0. 886 

b-6.406 

29/03/81 

5.49 

13/09/81 

6.20 

31/08/82 

5.34 

26/12/82 

5.60 

20/11/83 

5.36 

06/10/83 

6.03 

07/03/84 

5.52 

10/02/85 

5.69 

29/03/84 

6.00 

14/07/84 

6.18 

Table  II 
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Magnitude  computed  from  station's  factors 


station 

eO 

el 

e2 

S3 

S4 

e5 

e6 

e7 

e8 

e9 

LBF 

5.56 

6.15 

5.43 

5.54 

5.24 

6.03 

5.40 

5.88 

5.97 

6.17 

LOR 

5.50 

6.14 

/ 

5.62 

5.35 

6.03 

5.42 

5.92 

5.97 

6.17 

AVF 

5.66 

6.13 

5.17 

5.57 

5.37 

5.87 

5.61 

5.85 

6.02 

6.18 

LSF 

/ 

/ 

/ 

/ 

5.33 

/ 

5.48 

5.88 

6.01 

6.19 

SSF 

5.60 

6.12 

5.09 

5.63 

5.43 

5.91 

5.65 

5.90 

5.97 

6.12 

MFF 

/ 

6.15 

5.44 

5.50 

5.26 

6.02 

/ 

5.91 

5.98 

6.17 

GRR 

/ 

/ 

5.17 

5.63 

5.39 

5.92 

5.62 

/ 

5.97 

6.15 

FLN 

/ 

6.20 

5.31 

5.58 

5.33 

5.94 

5.47 

/ 

6.02 

6.17 

LDF 

/ 

/ 

5.24 

5.61 

5.37 

5.95 

5.51 

/ 

6.01 

6.17 

SMF 

5.63 

6.19 

5.39 

5.52 

5.28 

/ 

5.42 

5.87 

5.98 

6.17 

LRG 

5.67 

6.16 

5.20 

5.57 

5.33 

5.85 

5.59 

5.85 

6.02 

6.18 

TCF 

/ 

6.17 

5.25 

5.64 

5.29 

5.98 

5.55 

5.93 

5.98 

6.13 

RJF 

5.63 

6.18 

5.26 

5.54 

5.34 

5.92 

5.48 

5.87 

6.02 

6.18 

FRF 

5.63 

6.17 

5.30 

5.53 

5.31 

5.93 

5.47 

5.91 

5.99 

6.19 

LPO 

5.61 

6.17 

5.30 

5.53 

5.34 

5.95 

5.45 

5.87 

6.00 

6.19 

LFF 

5.58 

6.18 

5.26 

5.56 

5.34 

5.96 

5.48 

5.88 

5.99 

6.18 

EPF 

5.59 

6.16 

5.24 

5.57 

5.36 

6.05 

5.49 

5.82 

5.96 

6.17 

BSF 

5.60 

6.11 

5.14 

5.59 

5.38 

5.96 

5.58 

5.93 

6.00 

6.13 

HAU 

5.62 

6.16 

5.20 

5.58 

5.36 

/ 

5.50 

/ 

5.98 

6.16 

CDF 

5.55 

6.17 

5.15 

5.65 

5.37 

/ 

5.60 

5.90 

5.97 

6.10 

LMR 

5.64 

6.13 

5.22 

5.53 

5.31 

5.94 

5.57 

/ 

6.04 

6.15 

LPF 

/ 

/ 

5.23 

5.59 

5.37 

5.92 

5.56 

5.87 

6.02 

6.19 

CAF 

/ 

/ 

5.25 

5.62 

5.34 

5.99 

5.50 

/ 

5.99 

6.16 

MZF 

/ 

/ 

5.30 

5.53 

5.34 

5.93 

5.53 

5.89 

6.03 

6.19 

CVF 

/ 

/ 

/ 

5.57 

/ 

5.96 

5.54 

5.88 

6.01 

6.18 

BGF 

/ 

/ 

/ 

/ 

5.27 

/ 

5.55 

5.91 

5.98 

6.17 

HYF 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

5.90 

6.00 

6.18 

mean 

5.60 

6.16 

5.24 

5.58 

5.34 

5.95 

5.52 

5.89 

5.99 

6.17 

sig 

0.04 

0.02 

0.09 

0.04 

0.04 

0.05 

0.06 

0.02 

0.02 

0.02 

LDG 

5.5 

6.18 

5.3 

5.6 

5.3 

5.96 

5.5 

5.89 

6.02 

6.17 

Table  III 


Figures  Caption 


Fig  1  :  common  signals  for  the  ten  events  of  this  study 

Fig  2  :  Actual  signals  for  event  1  :  13/09/81 

Fig  3  :  Actual  signals  for  event  4  :  20/11/83 

Fig  4  :  Averaged  records  in  the  5  stations  of  Massif  Central 

Fig  5  :  Averaged  residual  signals  in  the  5  stations  of 

Massif  Central 
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The  set  of  seismic  data  is  the  one  already  used  in  our  previous  studies 
on  Lg  waves  {  18  quakes  with  ML  between  3*2  and  4.8,  within  the  crust; 
22  SP  stations  of  the  LDG  network  ). 

The  method  used  here  is  the  Algebric  Reconstruction  Technic  (ART) 
to  build  up  an  image  of  the  Q  factor  spatial  distribution. 

The  set  of  propagation  paths  is  presented  in  Figure  1. 

Each  point  of  the  .mage  of  Q»  takes  into  account  all  the  paths  within 
a  distance  of  50  km. 

In  order  to  apply  some  confidence  criterium  on  the  quality  of  this  Q 
evaluation  we  define  a  confidence  index  : 

i  c  =  n  source  x  1  rays 

where  : 

n  source  is  the  numbers  of  different  sources  which  generate 

the  rays 

1  rays  is  the  length  of  the  ray  inside  the  influence  circle  divided  by 
the  influence  radius. 

A  map  of  this  index  is  given  in  Figure  2,  where  the  represented  points 
fit  the  test  of  ray  length. 


-•'igure  2 


A*  a  re»ult,  all  the  blue  tonei  correspond  to  a  high  confidence  index, 
and  red  zone*  might  include  ayatematic  biaa. 

Figure  3  represents  a  rather  simple  regional  distribution.  If  one  juat 
keeps  within  high  confidence  limits,  the  main  regions  are  defined  as  : 

-  Limagne  Basin  ,  Forez,  Saone  valley  =  weak  figure  for  Q 

(  between  Q»  =  180  and  Q*  *  240  ) 

-  "Limousin  =  (between  Q-  =  200  and  Q»  =  260  ) 

-  South  Parisian  basin.  Berry,  Poitou,  North  of  Massif  Central 

(  between  Q°  s  300  and  Q  „  =  440). 

The  mean  Q»  value  (  always  computed  at  f  =  1  Hz  )  extracted 

from  these  results  is  :  Q<  =  290 

Similar  to  the  one  obtained  by  direct  Lg  measurement  (  Carr.pillo 
and  al  1984  ). 

We  need  to  extend  this  Q«  cartography  to  other  western  Europe  region 
by  using  other  stations  and  earthquakes. 


Conclusion 


These  results  pointed  out  the  importance  of  graben  zones  for  Lg  waves 
propagation.  Nevertheless,  at  this  stage  of  the  study  we  do  not  possess 
yet  any  separation  tool  between  two  types  of  phenomena  which  could 
create  attenuation 

-  small  scale  diffraction  (  or  scattering  ) 

-  large  scale  geometrical  effect 

A  study  of  frequency  dependancy  for  Q  should  help  us  in  the 
interpretation  of  our  data. 

On  the  other  hand,  developpements  of  numerical  simulation  methods 
encourage  us  to  investigate  more  complex  zones  in  order  to  interprete 
the  quasi  extinction  phenomena  of  Lg  waves  we  have  observed  along 
peculiar  paths. 


